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ABSTRACT 

We study the efficiency of growing large bodies, starting from a sea of equal-sized planetcsimal 
seeds. This is likely one of the earlier steps of planet formation and is related to the origin of the 
asteroid belt, the Kuiper belt and extra-solar debris disks. Here we study the case that the seeds do 
not collide frequently enough for dynamical cooling to be important (the collisionless limit), using a 
newly constructed conglomeration code, and by carefully comparing numerical results with analytical 
scalings. We find that large bodies grow primarily by accreting small seeds (and not by accreting each 
other). As the velocity dispersion of the small bodies (u) is increasingly excited by the growing big 
bodies, growth passes from the well-known run-away stage (when u is higher than the big bodies' hill 
velocity) to the newly discovered trans-hill stage (when u and big bodies both grow, but u remains 
at the big bodies' hill velocity). We find, concurring with analytical understandings developed in 
iLithwickl (|2013ft . as well as previous numerical studies, that a size spectrum dn/dR a JJ 4 results, 
and the formation efficiency, defined as mass fraction in bodies much greater than the seed sizes, is 
~ a few x Rq/cl, or ~ 10~ 3 at the distance of the Kuiper belt. This extreme inefficiency argues 
against the collisionless limit for the formation of the Kuiper belt and extra-solar debris disks. To 
conglomerate large bodies in those regions, the initial planetesimals may have to be small and highly 
collisional. 



1. INTRODUCTION 

Large solid bodies, ranging in size from less than a kilo- 
meter to thousands of km roam in the outskirts of both 
the solar and many extrasolar systems. In the case of the 
solar system, Pluto, Sedna a nd other large Kuiper belt 
objects are directly detected (iSlfpher fc Tombaughlll930l: 
Jewitt fc Lm3ll993l : [Brown et all 120051 : iSchlichting et al.1 
2009ft : and in the case of extrasolar ones, these bodies 



are inferred from the dust grains that are produced dur- 
ing their mutual collisions. Even a small amount of dust 
can intercep t enough star light to become detectable as 
debris disks (jAumann et al.lll984HWvattll2008| ). Surpris- 
ingly, ~ 20% of solar- type stars, even at an age similar 
to the Sun, harbor dust disks that are brigh ter than ours 
by m ore than three orders of magnitude (|Mever et al.l 
120071) . 

How do these large bodies form? Are they an interme- 
diate step in the formation route for planets? Why are 
they not incorporated into Earth or even Ncptunc-likc 
planets? Why are many exo-debris disks so bright yet 
ours so anemic? These are the questions we set out to 
answer. 

The conventio nal picture for the formation of our 
own Kuiper belt (IKenvon fc Luull99llOrmel et al.l201Cft 
ISchlichting fc Sarill201lL hereafter SS11) postulates that 
it started from ~ 10 Earth masses of primordial rock and 
ice (~ solid mass in the Minimum Mass Solar Nebulae), 
and only a fraction e ~ 10 -3 of this mass conglomerated 
into the large Kuiper belt bodies we see today. The rest 
was somehow removed. This extremely low efficiency oc- 
curs because, theorists argue, as bodies like Pluto grow, 
they stir the orderly orbital motion of the seed material 
to such an extent that the latter can no longer be accreted 



quickly into large bodies: Plutos starve themselves. 

However, there are four serious issues with this 
paradigm. First, exo-debris disks are too bright. The to- 
tal mass contained in their la rge objects must be ~ 1000 
times that in our own system (|Shannon !TWil20ll . Ac- 
cording to the conventional picture of Kuiper belt for- 
mation, this would have required ~ 10 4 Earth masses of 
primordial solids, some 1000 times more than the value 
expected in a minimum mass solar nebula. Second, a 
~ 10M e disk of kilometer bodies would have disrupted 
the long period Kuiper belt binarie s , which are observed 
in ab undance (|Parker et alj 120111 : iParker fc Kavelaarsl 
2012). Third, the removal of most of the primordial 
solids, by first grinding them down to dust grains and 
then blowing them out using radiation pressure, may 
not work on theoretical grounds (Lithwick et al. , in 
preparation). Fourth, the low formation efficiency is 
incompatible with some theories of Nept une migration 
(|Thommes et al.lll999t iGomes et all 120041) where a vast 
majority of the Pluto-like large bodies, if formed prior to 
Neptune's migration, were sti pulated to have been scat- 
tered away by giant planets (jSternl 119911 : iLevison et al.l 
12008ft . 

We explore a diff erent route, stimulated by the 
theoretical proposal of lGoldreich et a l. (|2004bL hereafter 
GLS) . Different from all p reviou s works (e.g.. iSafronov 
19691 iGreenberg etal] Il978t IWetherill fc Stewart 



19891 IKenvon fc Luul 119981: IKenvon fc Bromlevl 120081: 



Ormel et al. 1 120101 : ISchlichting fc Saril I2011D . we model 
the plantesimals as being so small that their mutual 
collisions are dynamically important. In this picture, as 
seed materials are stirred up their destructive collisions 
grind them down to smaller debris. Smaller seeds collide 
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much more frequently than larger ones and their relative 
velocities remain low. This collisional cooling guarantees 
efficient accretion by large bodies, until almost all mass 
has been converted to large bodies. Such a picture will 
allow the in situ formation of the cold classical Kuiper 
belt ([Shannon et al.ll2013a), and explain the hig h mass 
in bright exo-debris disks ( Shannon et al.ll2013bl ). 

But for the current publication, we restrict ourselves 
to the conventional case of no collisional cooling. There 
are three purposes to this paper. One is to document 
our coagulation code and certify its various components 
against analytical expectations (SJ3J). This code contains 
collisional break-down and cooling, but these two effects 
are unimportant here since wc adopt large initial sizes 
for the planctesimals. Second, our results demonstrate 
clearly (SJ3]) that, in the case of collisionless growth, the 
formation efficiency is limited to e ~ a few x a where 
a = -R©/a, i.e., a is the angle subtended by the Solar 
disk at distance a. In the Kuiper belt region, a ~ 1CT 4 . 
Hence only e ~ 10~ 3 of the initial mass could have been 
incorporated into large bodies that one observes there to- 
day (if growth was nearly collisionless) . Equivalently, to 
form these bodies the initial mass in planetesimals would 
have had to be ~ 10M®, similar to that expected in a 
MMSN disk extrapolated to the Ku iper belt distance, 
(|Weidenschillinejll977t lHavashil 1 1981ft ). Lastly, we com- 
pare results f rom o ur code against previous work (©. 
In particular, I S S 111 have recently proposed a simple ex- 
planation for the size distribution of Kuiper belt bodies, 
in the case of collisionless conglomeration. Wc examine 
their claim critically. 

This paper should be co nsidered a nume rical compan- 
ion to the recent work by ILithwickl (|2013|) where a con- 
glomeration phase called 'trans-hill growth' is first dis- 
cussed. That paper, largely analytical, arrives at similar 
conclusions as reached here. 

2. THE CONGLOMERATION CODE: COMPONENT 
TESTING 



We adopt a particle- in-a-box approach (|Safronovf l969) 
to study the interactions of plantesimals and their 
growth. The range of body sizes we track runs from 
~ km to - 10 3 kmQ Since the total number of parti- 
cles active in our simulations can exceed 10 15 , and we 
typically integrate the system for up to 10 6 dynami- 
cal times, we adopt by necessity a statistical approach 
where particles with similar properties (mass, eccentric- 
ity, inclination, semi-major axis) are grouped into the 
same bin. We do not yet have the capability of cou- 
pling N- body integrators with the s t atistical code, as was 
done b y iBromley fe Kenvonl (|2006f ) ; iKenvon fe Bromlevl 
(|2008D ; IGlaschkei ([20061 ) .As such, we fail to capture the 
evolution correctly at late times, when only a few large 
bodies dominate the evolution. 

2.1. Particles and Bins 

A total of six parameters are needed to specify an orbit. 
We choose the semi-major axis a, the eccentricity e, the 
inclination i, the argument of pcriapse w, the longitude 
of the ascending node O, and the mean anomaly M. To- 
gether with mass, each particle should be described by 

1 This is expanded in our next work when we consider grains as 
small as ~ fim. 



seven parameters. This is computationally prohibitive 
and approximations arc in order. In the following we ex- 
plain the simplifying assumptions we make in this code. 

We adopt a single- value of semi-major axis (a = 45AU) 
for all particles, spread to a width Aa ~ 0.13a ~ 6 
AU. Our single zone approach limits us to study sys- 
tems where the eccentricities are small (e < Aa/a). The 
particles are assumed to be distributed uniformly in w, 
tt, and M. We set i = 0.5e, as is approximately true 
if random energy is equipartitioned in each dimension 
(|Hornung et al] I1985D . So velocit y anisotropy, impor- 



tant when pa rticles are very cold (jlda fe Makinol [19921 : 
lRafikov|[2003D . is not properly treated here. 

These simplifications allow us to describe particles only 
by their masses and eccentricities. We track particle flow 
across bot h the mass bi n s and the eccentricity bins, as 
is done in IKrivoy et al.l (120051) . This differs from the 
approach of. e.g.. lKenvon fc Luul (|1998|) , where particles 
carry with them floating mass values, i.e., particles of a 
given mass have a single eccentricity. The eccentricity 
bins have values from e = 10 _7 toe=l and are equally 
spaced logarithmically with four bins per decade. The 
mass bins are also equi-spaced in logarithmic space, with 
six bins per decade of mass (i.e., Am w 1.47m). Particles 
of all masses are assigned the same bulk density, p = 
1.5g/ cm 3 . 

The dispersion velo city of a particle is 
(jWetherill fc Stewarti[l99l . 
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(1) 



where Wkep = ^o. with Q being the orbital angular fre- 
quency. 

Processes detailed in the following section drive mass 
and eccentricity evolution and we use a second-order 
Rungc-Kutta integrator to advance the simulation in 
time. There are two subtleties in the actual implemen- 
tation. First, when the outcome of an interaction re- 
sults in a new mass/eccentricity that falls in between 
two adjacent bins, we split the particles into the two 
bins conserving total mass and mass- weighted eccentric- 
ity (in the case of eccentricity), and total mass and total 
number (in the case of mass). Moreover, our integration 
has a fixed time step, typically set to be 100 years. In 
the case where dynamical evolution leads to undershoot- 
ing (negative number of particles) in a bin at the end 
of a timestep, we reduce the strength of the dynamical 
interaction such that the bin empties to zero. This is 
necessary when bodies are strongly affected by processes 
like viscous stirring and collisional cooling. 

2.2. Dynamical Friction & Viscous Stirring 

Mutual gravitational interactions lead towards 
equipartition of random kinetic e nergies between bodies . 
This is called dynamical friction ([Chandr asekhar 19431) . 
Additionally, in a keplerian disk, gravitational interac- 
tions between bodies converts their orbital energy into 
random kinetic energy. This is called viscous stirring 
(|Safronovlll969T) . 

To model dynamical f riction and viscous s tirring, we 
adopt the prescription of lOhtsuki et a l. ( 2002), who pro- 
vide semi-analytic formulae for the rates at which eccen- 
tricity and inclination evolve through gravitational scat- 
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Fig. 1. — : The absolute rates of viscous stirring (VS, 
red solid lines) and dynamical friction (DF, negative in 
sign and only acting on big bodies, blue dashed lines) 
used in our code (thick lines), as compared to the order- 
of-magnitude estimates (thin lines) of GLS, plotted here 
as functions of eccentricity (e). These rates are for a 
population of 1000 km big bodies interacting with a pop- 
ulation of small bodies, with equal surface densities, 
£ = a <~~j 0.1g/cm 3 , and equal dispersion velocities, 
u = v = eVkepj situated at 45 AU. The bend around 
e ~ 10~ 3 (Hill velocity of the 1000km bodies) occurs 
when the velocity dispersion transitions from sub-hill to 
super-hill. If one substitutes 2.5vh into the GLS expres- 
sions whenev er vh appears , as is suggested by N- body 
experiments (jNishidal 119831 : iGreenberg et al.l 119911 ) , the 
agreements between the two sets of curves will be much 
closer. 



terings, calibrated by N-body simulations. But when 
the mutual velocity of two bodies are below their mu- 
tual Hill velocity (sub- hill) , we ado pt instead the viscou s 
heating prescription in cq. (13) of iCollins et al.l (120071). 
calibr ated by the numerical simulations of lCollins &: Saril 

(pool . 

GLS give explicit formula for the rate of viscous stir- 
ring and dynamical friction, depending on whether the 
velocities are sub- or super-hill (also called shear- or 
dispersion-dominated). They consider two groups of par- 
ticles, smaller bodies with size s, surface density a and 
velocity dispersion u, and large bodies with radius R, sur- 
face density E and velocity dispersion v. The two most 
relevant expressions are viscous stirring of small bodies 
by the big ones, and dynamical friction of the big bodies 
by the small bodies. Assuming u > v, these are 




u > v H , 
u < Vh 

U > V H , 

u < vh 



(2) 



(3) 



Viscous stirring of large bodies by themselves can be 
found by substituting v for u in equation [2 As Fig. 
[T] shows, our adopted prescriptions agree with the above 
estimates well, when vh in the above expressions are sub- 



stituted by 2.5vh- Here, vh is the Hill velocity of large 
bodies, 



v H 



(M bie 
Ukep l3A/ ( 







1/3 



-1/3, 



SIR. 



(4) 



The big body mass Mbi g = 47r/3pi? 3 with p being the 
bulk density. We set p to be the mean density of the 
Sun, and use a to denote the angular size of the Sun, 
viewed by the bodies, a = Rq/cl. For the Kuiper belt, 
a ~ 10" 4 . 

2.3. Collisions: cooling and cascade 

Our conglomeration code includes the collisional pro- 
cess. Although not used in this collisionless study, it is 
essential for our follow-up studies of collisional growth. 
We document it here for completeness. 

In the particle-in-box method, the frequency of colli- 
sion for a particle is f c = nirb 2 v. Here, n is the number 
density of other particles, irb 2 the cross section for col- 
lision, and v the relative velocity. The relative velocity 
between two particles is simply set to be v = -\/T w i + wf) 
where Vi is as defined in equation ([TJ. 

Under our assumption of evenly distributed orbital an- 
gles, particles with eccentric ity e occupy a toru s in space 
that has a volume (see, e.g. lKrivov et al.ll2005l ). 



V(e) 



-in 



(l + e) 3 -(l-e) 3 



(5) 



Their number density is then obtained by n — N/V 
where N is the particle number. Between two groups 
of particles with different e's (ei,e2, but the same a), 
their overlapping volume is determined by the group with 
the smaller e (let it be ei), and the collision frequency 
is reduced by a ratio of V(ei)/V(e 2 ) to take account of 
the reduced residence time group 2 particles spend inside 
V(ei). 

We adopt a three-piece form (jGreenberg et al.l I1991L 
GLS) for the cross-section for collision between particles 
of size si and S2, 



7T (si + s 2 



1 



7Tb 2 = 



7T (Si 
7T (Si 



S2) 



S2 



V H < V, 
l 

a 2 vh < v < vh , 
i 

v < a 2 vh ■ 



(6) 

We follow GLS and name them super-hill, sub-hill, and 
super-thin cases. The last case occurs when the velocity 
dispersion is so small the particles can be regarded as 
infinitely thin in the vertical direction. These expressions 
apply when the velocity dispersion is isotropic (e ~ i). 

The outcome of collisions is modelled as changes in 
mass and eccentricity. For the eccentricity evolution, we 
set the orbits of the post-collision particles to be that 
for their center of mass, independent of the outcome of 
the collision. For head-on collisions, this corresponds to 
setting the coefficient of restitution to zero. Let the two 
bodies have initially masses mi, m 2 (m c = mi + 7722), 
orbital velocities v\, v 2 , the constraints of mass and 
momentum conservation lead to the follo wing expres- 
sions for the new orbital elements, (aj, e/) (|Krivov et al.l 
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Fig. 2.—: Collisional cooling, in our code (crosses) and 
from analytical estimate (solid curve). The eccentricity 
of particles are plotted as functions of rescaled time, tv, 
where v is the collisional frequency (eq. |HJ. Different 
symbols are for when s = 1 cm and when s = 1 mm. 
Numerical cooling rates agree with analytical estimates, 
as long as the latter rates are decreased by a factor of 5. 
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(7) 



where r is the position of the collision. In our single-zone 
calculation, we reset a/ to a after applying the above 
equations. In addition, we evaluate the vi • v 2 term by 
taking appropriate orbital averaging of the two orbits. 

An order-of-magnitudc estimate for the rate of colli- 
sional cooling among small bodies is (GLS) 



1 du 
u dt 



an 

ps 



(8) 



where v is the collisional frequency We confirm that our 
collisional code reproduces this scaling (Fig. [2]), albeit 
with a constant factor of ~ 5 offset. 

In terms of mass evolution, collisions can lead to catas- 
trophic destruction, (inelastic) rebound, or conglomera- 
tion. Catastrophic collision is defined where the primary 
body loses > 50% of its mass. We specify this to happen 
if the specific kinetic energy in the impact, 



roima „,2 



2 (m 1 +m 2 ) 

excee ds the disruption threshold ([Stewart fe Leinhardtl 
HH), 



Q* « 500 



-1/9^0.8 



+ 10- i (sl + S iy-v°-\ (9) 



with all numbers in cgs units. This is the scal- 
ing for weak aggregates, and may be appropriate for 
Kuiper belt bodies. For km-sized bodies, Q* ~ 



150 (v/1 



,0.8 



erg g 1 . This corresponds to a dis- 



ruption velocity of v ~ 100 cm/ s, or e « 10 -3 at 40 AU. 

When bodies of size si are catastrophically disrupted, 
we re-distribute their masses to smaller size bins with a 
number distribution that is power-law in size, dn/ds oc 
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Fig. 3. — : Evolution of the differential mass distribution 
during a collisional cascade, plotted at t = 0, 10 4 , 10 5 , 10 6 
yrs (red curves). Material strength is taken to be con- 
stant. As a result, one expects that bodies in collisional 
equilibrium satisfy dn/ds oc s -3 ' 5 (dotted line) In this 
simulation, initial surface density a = 0.1g/cm 2 , and 
mass of the catastrophically disrupted bodies are dis- 
tributed to smaller sizes as dn/ds oc s -4 . The final size 
spectrum is insensitive to this debris redistribution: the 
blue dashed curve shows the size distribution at 10 6 yrs 
if the redistribution instead follows dn/ds oc s~ 2 . We as- 
sume all particles smaller than l/_tm are instantly blown 
away by radiation pressure, which leads t o the wavy pat- 
tern n ear the cut-off size, as described in lThebault et al.l 
(|200l . 



s q for s < s±. We typically choose q 1 — —3.5. This 
distribution results in most of the mass ending up in the 
largest fragments. And Fig. [3] shows that results do not 
vary when we vary the value of q' from —2 to —4. Fur- 
thermore, to imitate the effect of radiative pressure, we 
remove all mass in size bins below our minimum size (at 
one micron). In the future, it may be useful to imple- 
ment the depen dence of the largest fragmen t mass on the 
impact energy (jStewart fc Leinhardt]|200"9l ). 

To verify the collisional mass evolution, we test our 
code against the standard case of collisional equilibrium 
for which an analytical solution is known. For a material 
strength Q* oc s p , the collisional cascade will carry a 
constant mass flux downward in partic le size and build 
up an equilibrium size distribution of (jDohnanvil 119691 : 
IDurda fc DermottllT997[ ). 



du 
ds 



oc s 



(10) 



with q = (21 + p)/(6 + p). The well-known Dohnanyi 
law is the special case with p = and hence q = 3.5. In 
the numerical experiment shown in Fig. [3l we adopt a 
constant strength of Q* = 10 8 crgs/g, and initialize the 
particle size distribution with dn/ds oc s . We ignore 
collisional cooling here by restricting bodies to only one 
eccentricity bin (e = 0.15). Irrespective of the power- law 
we choose for the debris redistribution, we find that small 
bodies in the system settle into the expected q = 3.5 
form and larger bodies gradually enter into collisional 
equilibrium as time goes on. 
If the collision energy is too low to cause destruction, 
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we separate the outcomes into two further categories: 
rebound or conglomeration. In our code, conglomeration 
occurs when the relative velocity falls below ten times the 
mutual escape velocity. This disfavors conglomeration of 
small bodies: kilometer-sized bodies only accrete each 
other when e < 1CP 3 at Kuiper belt distances. 

When this is not satisfied, the bodies rebound, with 
both having the same e (=e / , Eq. [7j) . We do not model 
either cratering collision or sticking by chemical forces. 
In the following, we discuss in more detail the accretion 
process. 

2.4. Collisions: accretion 

The collisional cross sections in eq. ([6]) lead to the fol- 
lowing rates of accretional growth (GLS), for bodies with 
size R and velocity dispersion u, accreting smaller bod- 
ies of surface density a and velocity dispersion u (with 

u > v) 



1 dR 
~R~~dt 



~pR ( 



, v H <u < v csc , 

[^f-) , a^v H < u < v H , 
i 

a 2 , u < a^vn ■ 



(11) 



We confirm that our numerical algorithm does reproduce 
the above relations. There is a subtlety here. These ex- 
pressions assume an isotropic velocity dispersion, i w e. 
However, when stirring is in the sub- hill regime (u < vh, 
eq. [2]), inclinations of bodies are excited at a lower 
rate than their eccentricities, while both quantities are 
damped by dynamical friction at comparable rates. We 
then expect these bodies to satisfy i <§; e and the accre- 
tion rate should proceed as the super-thin case (GLS), in 
which case the accretion rate is afl/(pR) ■ a -3 / 2 for all 
u < vh- In practice, however, stirring is contributed by 
bodies of various sizes, not just the largest sizes. Short of 
tracking the inclination evolution, it is difficult to ascer- 
tain the accretion geometry. So in this work, we simply 
assume an isotropic velocity dispersion even though it 
carries mistakes that may reduce accretion rate by up to 
a~ 1 / 2 (w/w/j). This, as it turns out, is more significant 
when dealing with accretion among big bodies. 

Here we comment on another numerical detail concern- 
ing the growth of the largest bodies. When two bodies 
(mi,m2, with mi < ra-i) merge, the combined mass in 
the new body typically falls between 7772 and the next 
mass bin (1713). We deposit two fractional bodies in the 
?7i2 and 7773 bins, respectively, requiring that the total 
mass be 777 1 + 777 2 and the total particle number be 1. 
This causes no major concern except where the fractional 
body in 7773 (number count 713) is the largest body in the 
system and it accretes much more quickly than other 
bodies. We then have a fractional body running away 
from the pack. To avoid this unphysical outcome, we 
adopt a so-called 'large body gating' procedure, which 
effectively allows the fractional body in 7773 to advance 
to the next mass bin (7774) only when the total mass in 
this bin has grown by 773 (7774 — 7773 ) . In such a gating pro- 
cedure, promotion occurs when the number of the largest 
body has reached at least unity. It emulates the growth 
speed of an individual large body, preventing inordi- 
nate run-away. Testing shows that thi s method repro- 
duces t he large body growth rates as in iGoldreich et al.l 
(|2004biT In reality, there is usually a single most mas- 



sive body in each accretion zone. Our approach here 
effectively assumes that the size of our simulation box 
is the size of the accretion zone. The validity of such a 
crude numerical approach should be verified with codes 
that incorporate N-body dynamics to reach beyond the 
statistical approach. 

3. CONGLOMERATION: RESULTS AND ANALYSIS 

In the remainder of this paper, we present the results of 
a number of collisionlcss simulations, in which rebound- 
ing and destructive collisions are disabled in the code, but 
conglomeration still occurs when bodies collide at a speed 
less than ten times their mutual escape speed. We focus 
on collisionless simulations in order to compare with pre- 
viously published results, in which typically small bodies 
collide less than a few times over the course of the sim- 
ulation. In an upcoming paper, we shall consider the 
effects of collisional damping and destruction. 

We present detailed results from a fiducial simulation 
in Fig. [U with Figs. [5] - |JJ] providing diagnostic details. 
For ease of comparison, we initalize o ur planetesimal 
disk s imilarly to those in previous works ([Kenvon fc Luul 
119981 SS11). Our disk is made up of s = 1 km planetesi- 
mals with a surface density of a = O.lg/ cm 2 , comparable 
to the solid value in a MMSN disk. Spread over a radial 
width of Aa/a = 0.13 at a = 45 AU, this corresponds to 
a total mass of IOMq. The bodies (~ 10 13 of them) are 
initially placed on dynamically cold orbits with e = 10~ 7 . 

The growth of the largest bodies is mainly due to ac- 
cretion of kilometer bodies, at all times. So an impor- 
tant quantity to focus on is eccentricity of the km bod- 
ies. The following scalings for the Kuiper belt region 
will prove useful: e csc = v CS c/wk C p ~ 10 _1 (i?/10 8 cm), 
where v csc is the escape speed from a body of size R, 
and e H = v H /v kcp ~ a 1/2 e osc ~ 10" 3 (i?/10 8 cm). We 
also discuss two related distributions: the number dis- 
tribution and the mass distribution. If the number dis- 
tribution is dn/dR cx R~ q , then the mass distribution is 
dm/dlogR oc i? 4 " 9 . 

The km bodies have a mean collision frequency 



v = 10 -7 yi- 
ps 



(12) 



This is the natural timescale in the problem of interest. 
So from now on, we express time always in unit of tv. 

3.1. Three Growth Stages 

The growth can be naturally divided into a few stages. 

1. Geometric accretion Initially, stirring is much 
faster than growth until e ~ e osc , and hence the initial 
growth is dominated by pair-wise conglomeration among 
equal-sized bodies (size s). Assume collisions occur with 
geometric cross-section. At time frCl, the average col- 
lision probability for each of the 10 13 km-sized bodies is 
{tv) , but the probability of a single body having experi- 
enced N collisions is (tv) N . The size distribution at this 
point falls off steeply with R. The largest body possible 
at time t is determined by (tv) N 10 13 = 1, or at a radius 
of 

-log(10 13 )i 1/3 



Ru 



log(^) 



•s. 



(13) 
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Fig. 4. — : Results from our fiducial simulation of collisionlcss conglomeration, starting from a disk of 1 km bodies 
with a surface density a = 0.1g/cm 2 , and an initial eccentricity of e = 10 -7 . The left panel shows the evolution of 
differential mass and the right panel eccentricity, plotted as functions of body sizes. Data are taken at 0.2, 4, 8, 12, 
16, 20 (all thin lines) and 40 Myrs (thick line), corresponding to tv = 0.02, 0.4, 0.8, 1.2, 1.6, 2.0 and 4. In the left panel, 
the dashed line indicates 1 body per mass bin, and the upturn in differential mass at the largest size bin is caused by 
the requirement that there be at least one full body in the most massive bin, as modelled by our 'large-body-gating' 
procedure (see W2A\ . At 40 Myrs, the number distribution of large bodies (100-1000 km, disregarding the upturn) can 
be roughly described by a power-law, dn/dR oc R~ q with q ~ 4. The efficiency of conglomeration, defined here as 
the mass fraction above 100km, is 10~ 3 . On the right panel, the two dashed lines correspond to the surface escape 
velocity and Hill velocity from bodies of that size. At 40 Myrs, small bodies are uniformly stirred by the largest bodies 
(e oc BP), while larger bodies experience dynamical friction from the seed planetesimals. 



So at tv ~ 0.3, the largest body has a size of R ~ 3 km 
(ignoring gravitational focusing; see below). By this 
time, km bodies have stirred each other to ~ 1/2 of their 
escape velocities (cq. [2]), or, e ~ 10~ 4 , much above en 
of the largest bodies already formed. 

2. Rapid Run-away If small bodies are super-hill, 
the largest bodies experience run-away. They pull away 
from the rest of the pack while carrying little total mass 
(left panel of Fig. [4]). Eq. (fTTj) indicates that the slowest 
growth occurs at the smallest size. So once the geometric 
accretion produces bodies with R ~ 3 x s, the bottle-neck 
is overcome and very large bodies can be conglomerated 
within a comparable amount of time. 

During this stage of rapid runaway, the small bodies 
dominate the stirring, and their eccentricities increase 
very gradually (sec Figure[5]). The big bodies grow faster 
than they can stir, and their size spectrum is very steep, 
q 4. With time, u/vh decreases towards unity and 
we enter the fol lowing newly-discovered phase of growth 
(ILithwickll20ll . 

3. Trans-hill growth (or slow run-away) Rapid 
run-away ends once the big bodies dominate the stirring. 
Thereafter, the big bodies' stirring interferes with their 
feeding, and the big bodies grow at the same rate as the 
small bodies' eccentricities are stirred. Equating the two 
yields 



vance criterion than what we adopt, for instance, the 
mass density S at a given size will be greater as the run- 
away (which thins out the number) is more restrained. 
So the transition would occur at a smaller R. This tran- 
sition point is also affected by the value of u. If the initial 
excitation is e = 10 -3 , for instance, the transition would 
occur at R = 10 3 km with again E/er = 10~ 4 . Lastly, this 
transition point is affected by the initial size spectrum of 
big bodies. 

Let us define a class of trans-hill bodies (-Rtrans (£)) 
which always satisfy vjj ~ u(t). These bodies are 
~ 100 km when the r apid run-away tr ansitions to slow 
run-away. Following iLithwickl (|2013| ). we argue that 
these trans-hill bodies remain effectively the main stir- 
rers in the system, J2 st i r ~ Rtrans, or equivalcntly, u ~ 
Uij(-Rstir)) in subsequent growth. Hence the name 'trans- 
hill' growth. 

To demonstra te this must be true (for a more detailed 
explanation, see iLith wick 20TT|), we consider the growth 
of bodies above and below -Rtrans (£)■ Growth of bodies 
with R < -Rtrans is in the super-hill (runaway) regime. 
Their growth rate (eq. Ill|) can be written as 



1 dR 

Rd{tv) 



Rt 



-a 1 x 



B 



B 



(15) 



S -1 (VH 
— rv 



1 



(14) 



For our simulation, this is satisfied when i? max ~ 100 km, 
with small body e ~ 10~ 4 , and E/cr ~ a ~ 10 -4 , at a 
time of tv ~ 1.2. 

The transition point is affected by our treatment of 
'large body gating'. If one employs a more stringent ad- 



Since this decreases with decreasing R, the size spectrum 
for R < i?trans is largely frozen. 

Growth of bodies with R > -Rtrans is in the sub- Hill (or- 
derly) regime. The shape of their size spectrum remains 
invariant with time. Since the initial size spectrum falls 
steeply beyond R > i?trans, stirring is initially dominated 
by bodies of size i?trans- And this will remain so through- 
out the trans-hill growth. 
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Fig. 5. — : Evolution of the eccentricity of km plantesi- 
mals is plotted here against the dimensionless time (solid 
red curve), for the simulation presented in Fig. |4] The 
Hill eccentricity of the largest bodies at any given time, 
eH(-Rmax), is plotted as a blue dashed curve In the geo- 
metric accretion phase (tv < 0.3), the small bodies self- 
stir to produce e oc i 1 / 4 . Rapid run-away follows once the 
growth bottle-neck is overcome. In this phase, e remains 
roughly constant as stirring is dominated by the ineffi- 
cient small bodies, while en rises steeply. After tv ~ 1.2, 
the growth of the big bodies is hampered by their own 
stirring, and the largest bodies remain trans-hill during 
growth, e ~ 2.5en After tv ~ a few, growth effectively 
stalls and e is again stirred up gradually as e oc t 1 / 4 . See 
text for more details. 



As a result of these two considerations, u sa t>H(^?stir) 
during the trans-hill growth. Equation (|14j) continues to 
hold, so the size spectrum extends to larger and larger 
sizes with E ~ const, or q rs 4. This is the value we 
fi nd in our simulationsn and is simi l ar to the results 
of lKenvon fcTuul (fllMUOrTnel et all (f20Toh . andlSUl 
An explanation for the slope of the size s pectrum, differ- 
ent from our above one, is given in ISSTTl . We discuss our 
different interpretations further in <2] 

In our simulations, bodies as large as 1000 km are 
formed by tv ~ 2. This is consistent with the mass- 
doubling time for the stirring bodies one obtains from 
eq. ([IT]), tv ~ aR/s. 

3.2. End of Growth 

The efficiency of formation, defined here as the mass 
fraction in bodies larger than 100 km, reaches ~ 10~ 3 
towards the end of trans-hill growth (Fig. [5]). The largest 
bodies formed are ~ 1000 km, or Pluto-sized. We argue 
below that the growth effectively stops when these values 
are reached. In other words, collisionlcss conglomeration 
produces large bodies but at a very low efficiency. 

During slow run-away, equation (|14p is continuously 
satisfied, yielding E/er ~ a. So the largest size bodies 
that can form in our simulation, insisting on at least one 
such body, is ~ 1000 km (see the line in Fig. [4] that 
represents each bin containing one single body). And 

2 The mass distribution shown in Fig. [4] shows a departure from 
5 = 4 for the largest bodies. This is a consequence of the fact that 
our numerical procedure insists on at least one single most massive 
body ('large body gating'). 
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Fig. 6. — : Fractions of initial mass that are incorpo- 
rated into large bodies (solid for bodies with R > 100 km 
and dashed for R > 10 km, as a function of i? max (also 
the most effective stirrer, i? s tir)- In the text, we define 
the efficiency of formation, e, as the mass fraction above 
100 km, with e ~ 10 -3 towards the end of the simulation. 



indeed we observe that evolution practically stalls when 
-ffmax ~ 1000 km. The formation efficiency, defined as 
the fractional mass accumulated in bodies greater than 
100 km, is then 



1 

e = — I Ed log R sa a few a 

a J 100 km 



(16) 



and is o f orde r 10 3 in the Kuiper belt region. In 
ILithwickl (I20T1 . we argue that trans-hill growth is typ- 
ically terminated when a single big body dominates its 
own feedi ng zone, i.e., when t he biggest bodies become 
oligarchs (jKokubo fc Idalll998l ). 

Our statistical code loses validity once the oligarchic 
phase commences. A statistical code coupled with N- 
body dynamics is needed. However, we can project the 
subsequent growth fairly robustly for the Kuiper belt re- 
gion, since oligarchy likely did not play an important role 
there. 

After trans-hill growth is terminated, small bodies are 
stirred to super-hill velocites, as u ~ t 1 / 4 until u reaches 
Vesc(-Rmax) (Fig- EJ- With such a hot population of small 
bodies, accretion by the big bodies proceeds with ge- 
ometrical cross-section. Yet geometric accretion is far 
too inefficient to form the observed population of Kuiper 
belt objects within the Solar system lifetime^ So there 
is no subsequent accretion to raise the efficiency of big 
body formation much above ~ 10~ 3 . This low forma- 
tion efficiency is a consequence of our neglect of collisions 
amongst small bodies. In a forthcoming paper, we show 
that when collisions are included, the formation efficiency 
can rise to nearly 100%. 

Growth is likely to evolve into the collisional regime. 
Even if plantesimals start as km-sized, as they are con- 
tinuously stirred, they approach their destruction speed 
at e ~ 10~ 3 (eq. O- Their mutual collisions beyond this 
point should be destructive and will produce many small 
debris. When this occurs, collisional cooling among small 

3 Pluto, for instance, has a geometrical optical depth of 
(1000km/40AU) 2 ~ 10" 1 ' 4 . 
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planctsimals can no longer be ignored. Growth will then 
be in the collisional limit. 

A minor caveat concerns our assumption of isotropic 
velocity dispersion for all bodies. For large bodies which 
stir each other in the sub-hill range, the inclination dis- 
persion may fall much below that of t heir eccentricities 
(|Rafikovi[200l iGoldreich et~aT]|2004bO . It is likely that 
accretion among large bodies could proceed then at a 
rate greater than that adopted here. 

3.3. Simulations with Different Initial Conditions 

Our above understanding of the conglomeration pro- 
cess allows us to explain the following dependences when 
initial conditions in the simulation are altered. In all 
cases, we find a similar size spectrum to that obtained 
abov e. We delay the discussion on the a-dependence to 

• initial eccentricity We compare the growth of 
bodies with initial eccentricities ranging from 10~ 7 
to 10~ 3 , with all other properties held constant. 
The escape velocity from km bodies is e ~ 10~ 4 at 
45 AU. We find that simulations with initial veloc- 
ity dispersion smaller than this quickly converge to 
our standard case; while those with velocity disper- 
sion above this value first undergo orderly growth 
and produce a more significant mass in intermedi- 
ate class bodies. In all case, the final efficiency is 
e « 1CT 3 . 

• initial surface density We vary the initial surface 
density of solids from 10 -3 times to 10 times the 
MMSN values (our standard case being 1 MMSN), 
with an initial eccentricity of e = 10~ 7 . We observe 
that the growth time scales inversely with the sur- 
face density (Fig. [7]). This is because the natural 
timcscalc in the problem is the small body colli- 
sio n time (eq. [T2| . This result has been reported 
bv lKenvon fc Luul (|1998f ). We also find that higher 
mass disks can harbor larger i? max - the final size 
of the largest body scales roughly as ct 1 / 3 , a result 
of us insisting that the largest bin has at least one 
whole body. Lastly, the efficiency of formation only 
depends logarithmically on the initial ct. 

• seed size We perform a suite of simulations with 
different starting sizes ranging from 10 m to 10 km. 
Since the mean collision time (eq. (fl2|) ) scales in- 
versely with starting size, and the entire evolution 
takes place within a few collision times, we expect 
that a smaller starting size leads to quicker growth. 
This is observed in Fig. [8] where the growth time 
scales roughly as v~ x ~ 10 7 (s/lkm) yrsQ The size 
of the largest bodies that form and the size spec- 
trum are comparable in all simulations. More im- 
portantly, the efficiency of formation (e) remains 
similarly low in all cases. 

To summarize, changes in surface density and seed size 
primarily affect the conglomeration results through their 
effect on the mean collision frequency (eq. [T2|) . 

4 This is s teepe r than the relation of s 1 /' 1 found by 
IKenvon & Luul Q199ST ). We suspect that this may relate to their 
initial conditions of super-escape velocity dispersions. 




Time (yrs) 

Fig. 7. — : Growth of i? ma x as a functions of time, for 
disks with different surface densities. The total mass 
spans from 10" 3 MMSN (bottom-most curve) to 10~ 2 , 
10~\ 1 and 10 MMSN (top-most curve), with the thick 
red curve being our fiducial case (Fig. |4]). We find the 
timescale of runaway growth obeys t cx ct -1 . All growth 
follow the pattern of a rapid run-away, followed by a slow 
run-away, and ending when the largest bin contains only 
one body. 




time (yrs) 



Fig. 8. — : Growth of the largest body in simulations with 
different seed sizes. Bodies grow slowly at first, with the 
velocity dispersion set by the escape velocity of the seed 
bodies. Eventually a few bodies run away, when they 
are sufficiently large. The timescale of runaway increases 
linearly with starting size. However, the largest bodies 
that formed are similar in size. 

4. COMPARISON WITH PREVIOUS WORKS 

4.1. IssTl 

While earlier work s (IKenvon fe Luul 119981 : 
IKenvon fc Bromlevl 120081 lOrmel et all I2010D have 
perform ed ex haustive studies of the conglomeration 
process, ISS11I is the first to present a simple analytical 
argument for the numerical results. Here, we critically 
e xamin e their argument. 

IssIH aims to explain the q = 4 spectrum, obtained by 
all works, as a result of equal accretion. They assume 
that the largest bodies grow by accreting each other and 
by accreting small bodies at equal rates ('equal accre- 
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tion'). This assumption, which was not properly justified 
or numerically confirmed in their study, enables them to 
derive the mass spectrum. They find <7 = 4, with a (con- 
stant) efficiency of large body formation of e ~ a 3 ' 4 . 
These claims appear to be backed up by their simple yet 
elegant numerical experiments. 

However, it is unclear why, on physical grounds, equal 
accretion should occur. Instead, we have shown that the 
q = 4 spectrum arises due to the fact that the largest 
bodies at any given time are trans-hill relative to the 
small bodies. We present a physical argument why this 
is naturally achieved ( £|3.1[) . As is shown explicitly in 
Fig. [9j we do not find accretion amon g large bodies to 
be important (also see iLithwickl 120131 ). Our expected 
efficiency of growth is ~ a (sec Fig. IT5j) . The difference 
between this and SSll's efficiency of a 3 / 4 is less easy to 
discern numerically since the two predictions differ by a 
very weak power of a. 
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Fig. 9. — : Cumulative contribution (by mass) to the 
accretional growth of i? max -sized bodies from bodies of 
various sizes, for the simulation shown in Fig. 21 when 
i? max = 30, 100, 300 and 1000 km. Accretion of small 
bodies always dominates the growth, with contribution 
from bodies of similar sizes reaching at best a few percent 
of the total accretion. This is in contrast to the 'equal 
accretion' assumption in lSSlll . 



For a closer comp arison , we have written a simplified 
code to emulate the ISS11I simulations. We adopt their 
set of initial parameters (4OM0 in 1km objects orbiting 
at 35 — 45AU with e = 3x 10~ 6 ), and their simplifica- 
tions whereby bodies feel dynamical friction from smaller 
bodies; bodies are viscously stirred only by bodies bigger 
than or equal to themselves; and all collisions lead to ac- 
cretion. We also adopt analytical dynamical friction and 
viscous stirring rates (equations [2l l3lll| l. For the collision 
rate, we adopt Equation (fTTjl Fl 

5 Our collision rate differs slightly from I S S 111 in the sub-hill 
regime, in that we ass ume acc reted bodies have isotropic veloc- 
ity dispersion (Eq. whereas ISS11I assume they are super-thin. 
This difference is unlikely to be important because in both our and 
their simulations the small bodies have u > vjj, while the biggest 

bodies are sufficiently dynamically cold (v < o^/ 2 vh) that the two 
rates are the s ame. Although we attempted to follow the exact 
prescription of I S S 111 , the sharp discontinuity in accretion rate at 



We compare results of our simulation to that of ISS11I 
in Fig. [10] and find good agreement in the formation ef- 
ficiency and the size spectrum. One difference, however, 
is that we find a drastically shorter timcscalc for growth 
in our si mulatio ns, by almost a factor of 100 compared 
to that in ISSlll . The results of this simplified simula- 
tion also qualitatively agree with those of our full coag- 
ulation simulation shown in Figure [4] — in particular, the 
size spectrum and efficiency at late times are compara- 
ble. However, there are a number of differences in the 
details, which are partly due t o the neglect of viscous 
stirring by smaller bodies in the ISSlll - type simulations. 

Focusing on the issue of equal accretion, Fig. HT1 shows 
that even in a simulation that is tuned to closely match 
ISSlll and that largely reproduces their results, accretion 
of small bodies dominates the growth. We do observe 
that accretion of big bodies stays at a constant fraction 
of the total growth (~ 10 %), as is the ca se for trans-hill 
growth (see Section 5.1 of lLithwickll2013l) . If we were to 
adopt the superthin accretion formula in our simulations, 
it is likely that this fraction would be closer to equal 
accretion. But the fact that our simulations produce q « 
4 shows that equal accretion is not a necessary ingredient 
for forming the size spectrum. 

Last ly, w e address the issue of formation efficiency. 
While ISSIH predicted that the final conglomeration effi- 
ciency should scale as a 3 / 4 , we predict e ~ a few x a. In 
Fig. [121 we explore a range of a from 10~ 2 to 10 -5 , corre- 
sponding to a heliocentric distance from 0.4 to 400 AU. 
We use the simplified code (rather than our full code) 
and find that our predicted scaling provides a good fit 
to th e numerical results, though we cannot exclude the 
ISSlll scaling due to the small numerical range. 

4.2. Kenyon & Luu (1998) 

We conduct a second comparison, to the simulations 
of iKenvon fc Luul (|1998f ). We choose this paper for com- 
parison over subsequent ones from the same group be- 
cause simulations in this paper experiences negligible in- 
fluence from the gas dynamics and have initial condi- 
tions that are the closest to our set-up. They begin with 
bodies of a single size (we will compare here to their 
800 m case), with a total mass of 10M® in an annulus 
from 32 to 38 AU and an eccentricity of e = 10~ 3 . This 
value corresponds to the escape velocity from 10 km-sized 
bodies. So the early accretion (before reaching 10 km) 
proceeds slowly without the benefit of gravitational fo- 
cussing. Adopting their initial parameters, our con- 
glomeration code yields size of the largest body at any 
given time, R ma , K {t), as in Figure ITB"! The growth time 
in our simulations ag ree within ~ 30% from those in 
IKenvon fc Luul (fl998l) . 

One significant differ ence appears w h en we compare 
the size distributions. IKenvon fc Luul (|1998l ) report a 
quantity r§ (defined by the cumulative number n(> r$) = 

v ~ vh for the superthin formula led to severe difficulties with the 
simulations. 

6 This difference in timescales is puzzling. But we note that 
when small bodies have u ~ vh, as they do in our simulations, 
the timescale to grow bodies of size R .is t ~ Rpct/(crQ) ~ 10 5 yr X 
(-R/100km) (Eq. Illll . whic h roug hly agrees with the timescale we 
sec in our simulation. In ISSllr s simulations, u ~ 3«h. which 
should lengthen the accretion time by the factor (u/vfj)^ ~ 10, 
not 100. 
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Fig. 10. — : Similar to Fig. [U but for a simulation that is tuned to closely match that of ISSlll . Results are plotted 
every 10 5 yrs, ending at 5 x 10 5 yrs. Both the number and size spectrum of large bodies agree roughly with those in 
ISSlll (blue dashed curve in the left panel). The timescales of growth, however, differ substantially. 
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Fig. 11. — : Same as Fig. HI but for a simp lified sim- 
ulation that is tuned to match that of ISSlll (see text). 
Instead of equal accretion of big and small bodies ex- 
pected by ISSlH we observe that accretion from small 
seeds dominates at all times, with accretion from large 
bodies staying at a small but constant fraction of the 
total growth, ~ 10%. 



10 5 ) of 51 km, while we find r§ = 31 km (roughly inde- 
pendent of time at late times). Correspondingly, our 
formation efficiency e is about 5 — 10 times lower. This 
difference may be explained by their pseudo-multizonc 
treatment where the big bodies' are prohibited from ac- 
creting each other when their respective Hill zones do 
not oy er lap. The same effect is observed in lOrmel et al.1 
(|2010f ) who also adopt multiple radial zones. We attempt 
to mimic this effect by imposing a minimum eccentricity 
floor of 10~ 4 to all bodies (our standard value is 10 -7 ). 
This keeps the big bodies artificially hot and prevents 
them from substantially accreting one another. We ob- 
serve a larger surface density in big bodies and we raise 
r,5 to 53 km in this case. 

5. CONCLUSIONS 




Fig. 12. — : The resulting differential mass spectrum from 
the simplified simulations, for different values of a, where 
a = 10~ 4 (40AU/a). The efficiency of formation e ~ 
a few x a (the horizontal lines denote 10a for different 
a value), as predicted by equation (TT51) . This is slightly 
steeper than the prediction of a 3 / 4 by ISSlll . 

In this contribution, we focus on the efficiency of form- 
ing large bodies by conglomeration, starting from a sea of 
small planetcsimals that are not cooled by frequent colli- 
sions. We find, concuring with previous studies, that the 
formation efficiency is e ~ a few x a, or ~ 10 ~ 3 at the 
distance of the Kuiper belt. 

We obtain this result numerically by constructing a 
conglomeration code that incorporates physical processes 
such as viscous stirring, dynamical friction and accretion. 
This code also has the capability of dealing with catas- 
trophic disruption and collisional damping. But in the 
current study with seed size of 1 km, collisional cooling 
is not relevant and we do not allow seeds to break down 
to smaller particles. 

We also deri ve the formation efficiency analytically 
(|Lithwickl 120131 ). We find that in collisionless environ- 
ments, growth passes through several stages, but natu- 
rally transitions to the most critical stage, which we term 
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Fig. 13. — : The size of the largest body as a function 
of time, usin g our conglomeration code (red curve) and 
from that of IKenvon fc Luul (|1998| ) (blue crosses). The 
growth rates agree qualitatively, although ours is faster 
by ~ 30%. 

'trans-hill' growth Growth is a run-away process when 
small bodies are super-hill and an orderly process when 
small bodies are sub-hill. Consequently, the biggest body 
at any given time sits at trans- hill, or u ~ UH-(i? max ). 
The fact that the system is driven to trans-hill growth is 
the key for why the size spectrum is dn/dR oc R~ q with 
q = 4, and why the formation efficiency is e ~ a few x a. 

To be confident of our numerical procedure, we have 
tested individual components in the code against the 
order-of-magnitudc formulae in GLS. We have also per- 
formed detail ed com par isons with previous sim ulations. 
These include ISSTfl and IKenvon fc Luul (fl998h . All pre- 
vious works yield a similarly low formation efficiency, 
and a similar size spectrum (q w 4) as we obtain here. 
But some detailed differences exi st. W e are able to re- 
produce the size distribution of ISS11I when we follow 
their numerical procedure, ignoring viscous stirring con- 
tributed by small bodies. However, we demonstrate that 
the key assumption they base their analytical under- 
standing on, that big bodies grow equally by accreting 
small bodies and by accreting other big bodies, is not 
what is responsible for the q = 4 spectrum and the effi- 
ciency. Large bodies grow primarily by accreting small 



seeds, and not by acc r eting each other. Comparison 
against IKenvon fc Luul (|1998l) show that we reproduce 
their growth timescale, albeit with a formation efficiency 
that is five to ten times lower. 

In the collisionlcss limit, trans-hill growth ends when 
the largest bodies become oligarchs, each responsible for 
stirring their own food. Subsequent accretion in this 
phase does not increase significantly the formation effi- 
ciency, at least in the Kuiper belt, due to the impossibly 
long accretion timescale for Pluto-type objects. 

Surface densities in large bodies in the Kuiper belt have 
been measured to be ~ 1CP 3 of that of the MMSN, and 
the largest bodies are ~ 1000 km in size. These have 
traditionally been viewed as successes for the collisionlcss 
coagulation theory. However, recent discoveries of bright 
extra-solar debris disks call this into question. Their dust 
luminosities reveal that they likely harbor large bodies 
that are a fa ctor of ~ 1000 in num ber than that in our 
Kuiper belt (jShannon fc Wull201lD . comparable to the 
total mass in a MMSN solid disk in these outer regions. 
This conflicts with the low formation efficiency, generic 
of collisionless conglomeration. 

In our simulations, km-size bodies reach such high ve- 
locity dispersion towards the end, that their mutual en- 
counter should cause fragme ntation into smaller p arti- 
cles. This is also observed in IKenvon fc Luul (|1999D . So 
even if planctesimals start as large as 1 km, collisional 
cooling may set in at some stage. An evolutionary path 
that is qualitatively different from that described here 
may ensue. 

To resolve the issu e of long forma t ion tim escale for 
Uranus and Neptune, iGoldreich et al.l (|2004a[) have pro- 
posed that conglomeration proceeds in a collisional envi- 
ronment, where small bodies are so small they are cooled 
by fre quent collisions. In future works (jShannon et al.l 
l2013af l . we follow this path and demonstrate that col- 
lisional conglomeration would also be able to raise the 
efficiency of formation to of order unity, thereby explain- 
ing, within one paradigm, the formation of the Kuiper 
belt and the extra-solar debris disks. 
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